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The  functional  differential  equation  Q' (t)  = AQ(t)  + 
BQi(T-t),  -°°  <t  < »,  whore  A,  B are  n * n constant 
matrices,  T _>  0,  Q(t)  is  a differentiable  n * n matrix 
and  QT(t)  is  its  transpose,  is  studied.  Existence, 
uniqueness  and  an  algebraic  representation  of  its  solutions 
is  given. 

This  equation,  of  considerable  interest  in  its  own  right, 
naturally  arises  in  the  construction  of  Liapunov  functionals 
of  difference  differential  equations  of  the  type  x(t)  = 

Cx(t)  + Dx(t-r),  where  C,D  are  constant  n x n matrices. 

The  role  played  by  the  matrix  Q(t)  is  analogous  to  the  one 
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TThe  functional  differential  equation  Q'(t)  = AQ (t)  + 

BQ  (x-t),  <t  < <»,  where  A,  B are  n * n constant 
matrices,  x >_  0 , Q(t)  is  a differentiable  n * n matrix 
and  QT ( t ) is  its  transpose,  is  studied.  Existence, 
uniqueness  and  an  algebraic  representation  of  its  solutions 
is  given. 

This  equation,  of  considerable  interest  in  its  own  right, 
naturally  arises  in  the  construction  of  Liapunov  functionals 
of  difference  differential  equations  of  the  type  x(t)  = 

Cx(t)  + Dx(t-x),  where  C,D  are  constant  n x n matrices. 
The  role  played  by  the  matrix  Q(t)  is  analogous  to  the  one 
played  by  a positive  definite  matrix  in  the  construction  of 
Liapunov  functions  for  ordinary  differential  equations. 

In  this  paper,  we  show  that,  in  spite  of  the  functional 
nature  of  this  equation,  the  linear  vector  space  of  its 
solutions  is  n2;  moreover,  we  give  a complete  algebraic 
characterization  of  its  solutions  and  indicate  computationally 
simple  methods  for  obtaining  these  solutions,  which  we  illus- 
trate through  an  example.  Finally,  we  briefly  indicate  how 
to  obtain  solutions  for  the  non-homogeneous  problem,  through 
the  usual  variation  of  constants  method. 
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played  by  the  matrix  Q(t)  is  analogous  to  the  one  played  by  a 

positive  definite  matrix  in  the  construction  of  Liapunov  functions 

for  ordinary  differential  equations. 

In  this  paper,  we  show  that,  in  spite  of  the  functional  nature 
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of  this  equation,  the  linear  vector  space  of  its  solutions  is  n ; 
moreover,  we  give  a complete  algebraic  character ization  of  its 
solutions  and  indicate  computationally  simple  methods  for  obtaining 
these  solutions,  which  we  illustrate  through  an  example.  Finally, 


ON  A FUNCTIONAL  EQUATION 


M 


ARISING  IN  THE  STABILITY  THEORY 
OF  DIFFERENCE-DIFFERENTIAL  EQUATIONS 

I . Introduction 

The  study  of  difference-differential  equations  has  received 
considerable  attention  in  recent  times  [2,6,7],  the  overwhelming 
interest  devoted  to  equations  with  positive  delays. 

In  this  brief  paper  we  wish  to  study  the  matrix  functional 
equation 

Q ' ( t ) = AQ(t)  + BQr(i-t),  < t < » (1.1) 

where  A,B  are  constant  n x n matrices  and  t > 0.  This  equation 

is  neither  of  the  retarded  nor  advanced  type.  We  show,  that, 

unlike  the  infinite  dimensionality  of  the  vector  space  of  solutions 

of  functional  differential  equations,  the  linear  vector  space  of 

2 

solutions  of  this  equation  is  n . Moreover,  we  give  a simple 
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algebraic  characterization  of  these  n linearly  independent 
solutions  which  parallels  the  one  for  ordinary  differential  equations, 
indicate  some  methods  of  computation  of  these  solutions  and  allude 
to  the  variation  of  constants  formula  for  the  nonhomogeneous 
problem. 

This  equation,  of  interest  in  its  own  right,  is  particularly 
important  since  it  arises  naturally  in  the  process  of  constructing 
Liapunov  functionals  for  retarded  differential  equations  of  the 
form  x’(t)  = Cx(t)  + Dx(t-t).  Datko  [4]  has  encountered,  in  a 


2 


somewhat  different  form,  this  equation,  but  has  not  studied  it. 

Repin  [11],  in  his  construction  of  Liapunov  functionals  uses  this 

equation,  but  does  so  erroneously  in  replacing,  in  (1.1),  the 
T 

term  BQ  (x-t)  by  BQ(x-t),  making  the  analysis  trivial. 

This  equation  has  been  used  by  Infante  and  Walker  [9]  in  the 
construction  of  the  Liapunov  functional  for  a scalar  difference- 
differential  equation.  The  study  presented  here  arose  in  the  use 
of  the  solutions  of  this  equation  in  a forthcoming  paper  [8] 
which  treats  the  construction  of  Liapunov  functionals  for  matrix 
difference-differential  equations . 

I I . Existence,  Uniqueness  and  Algebra ic  Structure 
of  the  Solutions 

Consider  the  equation 

Q'(t)  = AQ  (t)  + BQT(x-t),  < t < on,  (2.1) 

with  the  condition 

Q(j)  = K,  (2.2) 

where  K is  an  arbitrary  n x n matrix;  this  equation  is  in- 
timately related  to  the  differential  equation 

Q’(t)=AQ(t)+BR(t),  (2.3) 

1 < t < « , 

R'  (t)  = -Q(t)BT  - R(t)AT, 
with  the  initial  conditions 
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Q(j)  = K,  R[j)  = KT. 


(2.4) 


2 2 

Moreover,  for  any  two  n x n matrices  P,S,  let  the  n x n 
matrix  P 0 S denote  their  Kronccker  (or  direct)  product  [1,10] 
and  introduce  the  notation  for  the  n x n matrix 

s. 


s - (Sij)  = 


1* 


(s*x, • • • 'S*n^ ' 


■ th 


where  s . * and  s . . are,  respectively,  the  i row  and  the 

l * j 


.th 


column  of  S;  further,  let  there  correspond,  to  the  n x n 


2 T 

matrix  S,  the  n -vector  s = (s^ *,..., ) . 

With  this  notation  [1,10],  equations  (2.3)  and  (2.4)  can  be 

rewritten  as 


q (t) 

r (t) 

A 0 I 
-I  0 B 


B 0 I 
-I  0 A 


q(t) 
r (t) 


(2.5) 


and 


<j<i> 


= [k 


l*1 


r(2J 


T T T 

= ]kA^, . . . »k*nJ  , 


(2.6) 


which , with  the  obvious  correspondence  and  for  simplicity  of  notation, 
we  denote  as 


p (t)  = Cp  (t)  , (2.7) 

p(j)  = pT  . (2.8) 
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2 2 2 
Here,  p(t)  is  an  2n  -vector  and  C is  a 2n  * 2n  constant 

matrix. 

The  use  of  the  Kronecker  product,  which  has  allowed  us  to  re- 
duce (2.3)  - (2.4)  to  (2.7)  - (2.3)  permits  us  to  prove  our  first 
result . 

Theorem  1;  Equation  2.1  with  condition  (2.2)  has  a unique  solution 
Q ( t)  for  -ro  < t < 

Proof : If  equation  (2.1)  with  condition  (2.2)  has  a differ- 

T 

entiable  solution  0(t),  then,  defining  R(t)  = Q (t-r),  the 
pair  of  matrices  Q(t),  R(t)  will  satisfy  equations  (2.3)  and  (2.4); 
hence,  with  the  notation  introduced  above,  the  pair  of  vectors  q(t) 
and  r(t)  will  satisfy  equations  (2.5),  (2.6).  These  remarks,  the 

linearity  of  all  the  involved  equations,  and  the  uniqueness  of  the 
solutions  of  (2.5)  - (2.6)  immediately  implies  that  if  a solution 
Q(t)  exist,  it  is  unique. 

On  the  other  hand,  (2.5)  - (2.6)  has  a unique  solution  defined 
for  —on  < t r <»,  and  this  implies  the  existence  of  a unique  pair 
of  differentiable  matrices  Q(t)  and  R(t)  defined  for  -on  < t < ■" 
and  satisfying  (2.3)  - (2.4).  But  these  last  equations  can  be  re- 
written as 


^ Q(t)  = AQ  ( t)  + BR  ( t ) , 

d T T T 

R (x-t)  = AR  (T-t)  + BO  (T-t) , 


1 

(2.9) 


with  initial  condition 


Q(j)  = K = K (j)  , 


(2.10) 


from  which  it  follows,  from  uniqueness,  that  R(t)  = Q (r-t), 
completing  the  proof. 

Examination  of  the  above  proof  makes  it  clear  that  knowledge 

of  the  solution  of  (2.5)  - (2.6)  immediately  yields  the  solution  of 

(2.1)  - (2.2).  But  (2.5)  - (2.6)  is  a standard  initial  value 

problem  in  ordinary  differential  equations;  the  structure  of  the 

solutions  of  such  problems  is  well  known  [3,5].  Moreover,  since 
2 2 

the  2n  x 2n  matrix  C has  a very  special  structure,  it  should 
be  possible  to  recover  the  structure  of  the  solutions  of 
equation  2.1. 

Let  us  consider,  for  the  moment,  the  solutions  of  equation 

2 

(2.5).  Recall  [3,5]  that  it  has  2n  linearly  independent  solutions 

which  can  be  obtained  in  the  following  fashion.  Let  A^,...,A  , 

2 

p = 2n  , be  the  distinct  eigenvalues  of  the  matrix  C,  that  is, 
solutions  of  the  determinental  equation 


det [ AI-C]  = det 


( \ I -A ) 0 I 


I 0 B 


■B  0 I 


I 0 (AI+A) 


(2.11) 


each  A.,  j = l,...,p  with  algebraic  multiplicity  m.  and  geometric 

3 r s r 2 2 

multiplicities  n . , y n.  = m.  , ) m.  = 2n  . Then,  2n  linearly 
1 r=l  3 J j J 


independent  solutions  of  (2.5)  (or  (2.7))  are  given  by 


(t)  = e 

3 ,r 


A . (t  - i)  q (t  - ~ 


r.  q-i 


l 

i=l 


2 


(q-i)  ! j / r 


(2 


where  q = r = l,...,s,  l \ m_.  = 2n  , and 


r=l  ^ ^ • 3 


the  2n  linearly  independent  eigenvectors  and  generalized  eigen- 
vectors are  given  by 


[A . I-C] e 

3 l.r 


i-1 

-e  . , 

Jtr 


0 


e . = 0 , 

: » s 


(2, 


A change  of  notation,  and  a return  from  the  vector  to  the  matrix 
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form,  shows  that  2n  linearly  independent  solutions  of  (2.3) 
are  given  by 


p? 

3,r 

(t) 

(t  - 1)^1 * 3 

L1 

3 

\J/*3 

j rr 

(t) 

i=l 

(q-i) ! 

3 / r 
__  __ 

(2, 


for  q—  1, . . . , n . , r l,...,s,  j — l,...,p,  ) n . m. 

3 n=l  3 3 


r 2 i i 

) m . = 2n  , and  where  the  generalized  eiqenmatrix  pair  (L.  ,M. 

k 3 ^ 1 J,r'  j, 


associated  with  the  eigenvalue  A . , satisfy  the  equations 


(A . I-A) L1  - BM1 
3 3>*  j,r 


-Li_1 

3'? 


i T i T 

L B + M ( A . I+A  ) = 

3 i r 3 / r 3 


-M1-1  , 

3 ' r 


1 


.12) 


13) 


14) 


) 


15) 


w 


c , r ..  0 „0  . 

for  x = r = L.  = M.  = 0. 

3 3 » s i , s 


The  structure  of  these  equations  is  a most  particular  one; 


indeed,  if  they  are  multiplied  bv  -1,  transposed,  and  written  in 


reverse  order,  they  yield 


.T  .T  . T 

(- A • I -A)  M 1 - BL1  = M1  , 

1 ]»r  3,r  ],r 


(2.16) 


. T ,T  .T  . '1 

111  T 1 — 1 

M.  B + L (-A .I+A1)  = L. 

]-r  3 ],r 


T T 

r 0 0 

for  l = l,...,n.,  r = l,...,s,  L.  = M.  = 0.  But  this  result 
3 3 ,r  ] , r 


demonstrates  that  if  A ^ is  a solution  of  (2.11),  -A,.  will  also 


be  a solution;  moreover  A.  and  -A.  have  the  same  qeometric 

3 3 


multiplicities  and  the  same  algebraic  multiplicity.  Hence,  the 


distinct  eigenvalues  always  appear  in  pairs  (A.,-A_.),  and  an 


examination  of  equations  (2.15)  and  (2.16)  shows  that  if  the 


generalized  eigenmatrix  pairs  corresponding  to  A . are  (L.  ,M-?  ) , 

3 3 ' r 3 r £ 


the  generalized  eigenmatrix  pairs  corresponding  to  -A^  will  be 
((-i)1+1m.J-T  , (-1)  1+1liT  ) . 

J ir  3 ’ 


But  these  remarks  imply  that  if  the  solution  (2.12)  correspond- 


ing to  A.  is  added  to  the  solution  (2.12)  corresponding  to  -A^ 


multiplied  by  (-1)C^+'*',  the  n2  linearly  independent  solutions  of 


(2.3)  given  by 


B 


“3rS  X.(t  - J)  q (t  - j)q  1 Lj,r 

q C ili  Tq-TT! 

n?  _(t)  1 1 m1 

J/S  j , r 


-A . (t  - 


Tv  , . T.q-X 

? v ^ ~ 2 q+  i 

) r— 4-rn (-l)q 


i=l  (q"i)! 


L1 

1 . r 


qT 

t - n . i — 


satisfy  the  condition  r ( j)  = n j ^ r ('2')  - But  this  is  precisely 


condition  (2.4);  it  therefore  follows  that 


:q  (t)  = J-T [e  j (t  + (-l)q+ie  Aj(t  2>M?-,r],  (2.17) 

1 ,r  (q-i)  ! l/i  J/1 


i-i  (q-1,! 


for  q = l,...,nr,  f nr  = m . , J m.  = 2n  , are  n linearly  in- 
1 r=1  1 3 2j  J 

dependent  solutions  of  (2.1).  We  have,  therefore  proven 


Theorem  2:  Equation  2.1  has  n linearly  independent  solutions  given 

bv  equation  2.17,  where  the  generalized  eigenmatrix  pairs  (L^  ,M*  ) 

satisfy  equation  2.17  for  one  of  the  elements  of  the  pair  (Ay-A^), 
each  of  which  is  a solution  of  equation  2.11. 


It  is  interesting  to  remark  that  the  determinental  equation 

2 2 . 
(2.11),  involving  a 2n  x 2n  determinant,  given  the  commutativity 

of  its  elements,  can  always  be  rewritten  as 


det [ ( A I-A)  0 ( Al+A)  + D 0 D]  = 0. 
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Theorem  2 gives  the  desired  algebraic  representation  of  the  solu- 
tions, a representation  which  is  completely  analogous  to  that  for 

ordinary  differential  equations.  It  is  surprising  that  the  vector 

2 

space  of  solutions  of  (2.1)  has  dimension  n . 


Ill . Some  Further  Characterizations 

Theorem  2 of  the  previous  section  gives  a complete  characteriza- 
tion of  the  solutions  of  our  original  functional  equation.  It  is 
possible,  however,  to  give  some  further  properties  of  the  eigen- 
matrix  pairs  in  certain  particular  cases;  these  further  characteriza- 
tions of  the  eigenmatrix  pairs  are  very  useful  from  a computational 
viewpoint,  as  we  demonstrate  in  the  next  section. 


Lemma  1 : If,  for  a X^  satisfying  equation 

an  Uj  / 0 and  corresponding  vectors  xj'^j 


2.18,  there  exists 
such  that 


[ X . I - A - a . D ] x . = 0 , 

1 1 1 

[X  .1  + A + — B ] y . = 0, 
J 3 


(3.1) 


1 T 

then  equation  2.15  is  satisfied  for  i = 1 with  L.  = x.y.  , 

J,r  JrJ  ' 

n1  = m.L1 
m , r j j , r 


Proof : A simple  substitution  of  the  indicated  result  into 

equation  (2.15)  shows  that  this  equation  is  satisfied. 

The  proofs  of  the  following  two  lemmas  are  equally  obvious. 
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Lemma  2;  If,  for  a Aj  satisfying  equation  2.18,  there  exist 
vectors  xj'^j  such  that 


Bx  . = 0 
3 


(3.2) 


(A . I+A) y j = 0 


then  equation  (2.15)  is  satisfied  for  i = 

,,1  T 

M . = x . y . . 

D/r  r: 


1 with  L = 0, 

1 r r 


Lemma  3 ; If,  for  a A satisfying  Equation  2.18  there  exist 
vectors  Xj,yj  such  that 


(A . I -A) x . = 0 
1 1 


(3.3) 


By . = 0 
' J 


then  equation  2.15  is  satisfied  for  i = 1 with 


1,1  = x . y1  , = 0. 

I ,r  \J  l ' ) , r 


It  is  noted  that,  if  the  assumptions  of  any  of  these  three  lemmas 
hold,  then  the  form  of  our  eigenmatrix  pairs  is  dyadic.  More- 
over, determination  of  the  depends  on  simultaneous  solutions 

of  the  determinental  equations,  for  a given  A ^ satisfying  (2.18), 


dot  I A j I-A-u ■ U ) = 0, 


dot [ ' . T + A + — Hi  = 0, 

"j 


(3.4) 


L. 
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where,  for  Lemma  2 we  let 


1 


a . 


1 


0 and,  for  Lemma  3,  u. 


0.  The 


computations  involved  for  the  determination  of  a.,x.  and  y are 

3 3 3 

much  simpler  than  those  implied  by  equation  (2.15). 

If  B and  either  (A^I-A)  or  (A^I  + A)  are  not  invertible, 
then  Lemmas  2 and  3 apply.  it  is  not  difficult  to  show  that,  if 
such  is  not  the  case,  there  always  exists  at  least  one  j-  0 

such  that  (3.4)  is  satisfied,  which  implies  the  applicability 
of  Lemma  1. 


Lemma  4:  If  either  B is  invertible  or  (A.I-A) 

1 

are  invertible,  then  there  exist  an  a . / 0 which 


and  (AjI+A) 
satisfies  (2.22) . 


Proof : If  B is  invertible,  for  a given  j,r  and  i = 1, 

equations  (2.15)  are  equivalent  to 


- 1 1 T T ± 1 

B (A.I-A)L.  (A.I+A)  B = -L 

3 J,r  j j,r 

(3.5) 

-1  ] T rp~  1 1 

B (A  . 1 -A)  M . (A  . HA)  B = -M  . 

3 J » r j j , r 


But  this  implies  that  the  matrix  [B  1 ( A I —A ) ] O [B  1(A..I+A)] 
has  at  least  one  eigenvalue  that  equals  -1.  Since  [10]  the  eigen- 
values of  a Kronecker  product  are  the  products  of  the  eigenvalues  of 
the  two  matrices  appearing  in  the  product,  then  there  exist  an 

a • / 0 and  vectors  x . and  v • such  that 
3 3 3 


-1 


(A.I-A)x.  = a X 
J J .1 


1 


(3.6) 


B ( >.  . J+A)  v . 
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which  is  equivalent  to  the  assumptions  of  Lemma  1. 

If  B is  not  invertible  but  (A^I-A)  and  (A^I+A)  are  both 
invertible,  then  (2.15),  for  i = 1,  are  equivalent  to 


( A . I-A) _1B  L.  BT(A.I+AT)_1 

> ) r L J 


-L 


1 

j,r 


- 1 1 T T - 1 

(A . I-A)  B M7  B (A . I+A  ) 

D ) i r 1 


-M1 


1 


and  a repetition  of  the  above  argument  leads  to  the  same  conclusion. 

These  last  four  Lemmas  imply  that,  associated  with  each  distinct 
eigenvalue  pair  (Aj,-Aj),  a solution  (2.17)  to  our  equation  exists 

in  which  If.  and  M'.  are  for  some  r dvadic  and  linearlv  dependent. 
It  is  also  worth  remarking  that  Lemmas  1-3  are  not  exclusive,  and 
that  more  than  one  set  could,  in  some  cases,  be  ob- 

tained from  the  appropriate  equations,  yielding  a similar  result 
for  other  values  of  r. 

IV.  An  Example 

Here  we  present,  as  an  illustration  of  the  computations  in- 
volved, the  construction  of  the  linearly  independent  solutions  for 
a 2 x 2 systems  of  the  form 


T 


Q'(t)  = AQ(t)  + BCV(x-t),  — • < t < ••  . 


(4.  i; 


v/here,  for  illustrative  purposes,  we  have  chosen 


A = 


-5  0 


B = 


0 -2 


D 


4 0 

2 0 


(4.2) 


For  the  computation  of  the  eigenvalues,  the  determinental  equation 
(2.18)  gives 


do 


( X I - A ) 0 1 -B  0 I 


I 0 B 1 0 01  + A) 


det|(\I-A)  0 (M  + A) + B 0 B]  = 


- det 


(A+5) ( A-5 ) + 16  0 0 

8 (X  + 5)  ( A — 2 ) 0 

8 0 ( A + 2 ) (A-5) 

4 0 0 


0 

0 

0 

(A  + 2)  ( A-2 ) 


= 0, 


from  which  we  obtain  the  four  pairs  of  roots  (A^,-A^)  = (5,-5), 
(A2,-A2)  = (3,-3),  ( A 3 , - A 3 ) = (2,-2),  (A4,-A4)  = (2,-2)  = (A3,-A3). 

For  the  first  pair,  note  that  A^  = 5 is  an  eigenvalue  of  -A 
and  that  B is  not  invertible.  Application  of  Lemma  2,  equation 
(3.2)  yields  that 


Bx  = 0,  or 


4 0 

2 0 


x3  = 0, 


( A^ I+A) = 0,  or 


0 0 
0 3 


= 0 ; 


hence  we  have  L j ^ = 0 and  m|  ^ 


1 [1,0]  yielding  the 


sol ut  ion 


- -- 


The  second  pair,  (>,.,, -\0)  = (3,-3),  does  not  consist  of  eigen- 
values of  either  A or  -A.  Hence,  we  search  for  an  a such 
that  equations  (3.4)  are  satisfied,  that  is 


which  yields  the  value  ~ 2.  Equations  (3.1)  with  X j — 3 and 

u . = 2 now  yield 
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i 


The  last  pair  of  eigenvalues  is  of  algebraic  multiplicity 
two  but  it  is  easy  to  check  that  they  are  of  geometric  multiplicity 
one,  hence  we  can  attempt  to  treat  them  once  again  through  our 
lemmas . 

In  the  first  case,  '3  =2  is  an  eigenvalue  of  -A.  Then 
equations  (3.2)  become 


Bx  3 1 = 0 ' or 


4 0 

">  2 


*3,1  = °' 


(\3I+A)y3  l = 0,  or 


-3  0 


0 0 


y3,l  = °' 


yielding  L 3 ^ = 0 and  ^ = 

solution 


[0,1],  and  therefore,  the 


. 1 

“3,1 


(t) 


e"2(t-s) 


0 

0 


0 

1 


(4.5) 


In  the  second  case,  V , - +2  is  an  eigenvalue  of  A,  still, 
equations  (3.4)  are  satisfied  for  an  ^ 0;  indeed,  we  require 
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det[\  I + A + — B]  = det 
3 *3 


2 - 5 + — 0 


2-2 


which  are  clearlv  satisfied  for 


u 3 4 * 


The  corresponding  vectors  are  immediately  computed  as 


X3,2 


yielding,  through  Lemma  1,  the  solution 


1 (t)  = 2 (t-3) 

“3,  2 11  e 


0 

8 

+ I 

o 

0 

0 

7 

4 

8 

7 

-2  (t-3) 


It  is  easily  seen  that  these  four  solutions  are  linearly  in- 
dependent. 


V . The  Honhomogenoous  Problem 

We  briefly  consider  the  nonhomogeneous  problem 


Q' (t)  = AQ ( t ) + BO  (x-t)  + F ( t ) , < t 


where  A,B  are  constant  n x n matrices  and 


F(t)  is  a continuous 
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n x n matrix.  We  seek  a particular  solution  of  this  problem. 

2 

For  simplicity  of  notation,  let  the  n linearly  independent 

solutions  given  by  (2.17)  be  relabeled  as  Z,  (t)  = (z,  • (t) ) , 

K -L  / J 

0 

k = l,...,n".  Then,  in  a manner  completely  analogous  to  that  for 
ordinary  differential  equations,  we  obtain 

Theorem  3 : Equation  5.1  has  a particular  solution  given  by 


n 


2, 


q (t)  = y 

k=l 


r^ (s) ds 


zk( 


t)  , 


(5.2) 


where  the  r^(s)  are  cJiven  below. 


Proof : Substitution  of  (5.2)  into  (5.1),  given  that  the  Z (t)  are 

solutions  of  the  homogeneous  equation,  yields  that  the  r^(s)  must 
satisfy  the  equation 


)'  r.  U)zk(t)  - F ( t ) . 
k-  1 K 


(5.3) 


In  the  notation  of  Section  2,  define  the  n -vectors 


, k 


k 


T 


r j.  ( t ) = (z(t)^,  . . . (t)  ) , r (t)  = ( r x ( t ),...,  r 2 ( t ) , 

f (t)  = (f  f *)',  and  the  n“  x n“  matrix  ®(t)  = 

1*  n* 

( r , . . . , r ).  Then  (5.3)  is  equivalent  to 
n*" 


(h)  ( t ) r ( t ) = f (t)  . 


(m)(t)  is  clearly  nonvanishing  for  all  t,  given  the  linear  independence 


of  the  solutions  Z'(t);  hence 


i 


- 


13 


r(t)  = (H)-1(t)f(t) 
and  this  concludes  the  proof. 


■■ 
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